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Abstract. - Using Monte Carlo simulations and a mean-field theory, we study domain growth 
in a driven lattice gas. Mediated by a single vacancy, two species of particles ("charges") 
unmix, so that a disordered initial configuration develops charge-segregated domains which 
grow logarithmically slowly. An order parameter is defined and shown to satisfy dynamic 
scaling. Its scaling form is computed analytically, in excellent agreement with simulations. 

Introduction. Driven diffusive lattice gases, first introduced by Katz, et.al. display 
generic non-equilibrium features in their whole phase space Even when interactions are 
due to excluded volume constraints alone, complex phase diagrams can be induced by local 
dynamic rules ||-||] which violate detailed balance Q| . While disordered phases and universal 
properties near continuous transitions are quite well understood by now, ordered phases have 
proven far more complex In particular, studies of phase ordering and coarsening in driven 
systems have revealed surprising behaviors, quite distinct from those 0| exhibited by systems 
evolving towards an equilibrium final state. For example, when driven, the two-dimensional 
(2D) Ising lattice gas develops non-universal domain morphologies which grow in a highly 
anisotropic fashion and do not satisfy dynamic scaling [p|. I n ID, ordered domains grow with 
time as t^^^ |P, in contrast to the equilibrium t^^^ law |lO|. Analytic results are sparse and 
focus primarily on establishing growth laws [|[ ^, |l2| . Clearly, further studies of domain 
growth far from equilibrium are needed before a more general framework can be established. 

In this letter, we report a detailed dynamic scaling analysis for defect- mediated domain 
growth in a driven three-state lattice gas (|]. Two species of particles, labelled "positive" 
and "negative", diffuse on a periodic lattice by hopping onto a single empty site (the defect). 
Only nearest-neighbor jumps are allowed, biased by an "electric" field E aligned with a lattice 
axis. At finite vacancy concentration, this system exhibits a phase transition, controlled by 
drive and particle density, from a disordered, "free-flowing" phase to an ordered, "jammed" 
phase Ordered configurations are spatially inhomogeneous: positive and negative 

particles form adjacent strips transverse to E, impeding each other. A mean-field theory 
predicts the structure and stability limits of both phases Models of this type have 

been invoked to describe water-in-oil microemulsions Jl6t , gel electrophoresis JT?! and traffic 
flow @. 

Remarkably, even the presence of a single vacancy suffices to order an initially random 
system. This problem, in both its static and dynamic aspects, is an example of a much-wider 
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Fig. 1 - Snapshots of a 16 x 24 lattice at E = 0.8, taken after (a) 10^, (b) lO", 
lO'^ MC steps. Negative (positive) particles are colored black (white), and the hole is gray. E points 
upwards. The initial configuration was random. 

ranging class of interacting random walk and defect-mediated domain growth problems. The 
vacancy plays the role of a highly mobile defect [ p^ , whose motion restructures its environment 
(i.e., the particle configuration), while the latter provides a feedback through the local .jump 
rates. In the resulting steady state, investigated in detail in particles form a charge- 
segregated strip, with the vacancy localized at one of the interfaces. The density profiles obey 
characteristic scaling forms, controlled by drive and system size. 

Here, we focus on the time evolution of this system in 2D, seeking to understand how 
ordered domains form, grow and finally saturate. We present Monte Carlo (MC) data, sup- 
ported by an approximate solution of the time-dependent mean-field equations. The ordering 
process exhibits three regimes: the initial formation of the strip, an extended logarithmic 
growth regime and the crossover to steady state. We demonstrate that the system satisfies 
dynamic scaling in the second regime and compute the scaling function analytically, in excel- 
lent agreement with the data. We conclude with some suggestions for further study. More 
details can be found in |^ . 

The microscopic model. Our model is defined on a 2D square lattice of x Ly sites with 
fully periodic boundary conditions (PBC). The occupation of each site is represented by a 
spin variable Sxy taking the values +1, —1,0, if a positive or negative charge or the vacancy 
is present at site {x,y). For simplicity, we restrict our study to equal densities of positive and 
negative charges. At each MC step (MCS), the vacancy exchanges with a randomly chosen 
nearest neighbor Sxy, with Metropolis rate |^^, min {1, exp{sxyE6y)}. Here, Sy = 0, ±1 is the 
change of the particle's j/-coordinate due to the jump. E denotes the bias, miiform in space 
and time and directed along positive y. 

This deceptively simple dynamics induces a charge segregation process on the lattice, as 
illustrated by Fig. 1. Starting from a random configuration (not shown), the system remains 
disordered for early times (Fig. la). Eventually, the hole begins to segregate the two species: 
In Fig. lb, an interface between regions of opposite charge begins to develop. Due to PBC, 
a second interface must also form; this occurs at a time set by Ly. After 10^ MCS (Fig. Ic), 
the segregation of charges is already quite apparent. Clearly, the field drives the hole towards 
the lower (the "downstream"), and away from the upper (the "upstream"), interface. Typical 
configurations are homogeneous in the transverse direction. For our choice of parameters, 
the interfaces are well separated [Q. Since any increase of order requires a series of field- 
suppressed exchanges, the approach to steady state is very slow (Figs. Ic-e). 

In the following, we provide a brief quantitative analysis of this process, beginning with 
our MC results. The control parameters of our study are E, ranging from 0.2 to 2.0, and 
the system size, = 16 or 40, 16 < Ly < 72. The initial configuration is random. Time is 
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Fig. 2 - Charge density profiles for Lx = Id, E = 0.8, and a range of Ly, at different MC times 
(arrows). Some data points, away from tfie central region, are not sfiown for clarity. 

measured in MCS. All of our data are averaged (denoted by (•)) over 100 independent runs 
(samples), resulting in good statistics with errors well below 5%. 

To probe the growth of the ordered domain, we measure averaged hole and charge density 
profiles, defined via, (l){y,t) = (^Y~J2xi^~^xy)'^ and^(2/,t) = (^-^J2x^xy^- Due to trans- 
lational invariance, strips can be centered at any y; thus, care must be taken when averaging: 
In each sample, we determine the maximum of the hole density over a sufficient time interval. 
This maximum marks the downstream interface. The charge density profiles from different 
samples are now shifted so that these maxima coincide, and averages can be taken. Clearly, 
this procedure is not very reliable for early times, since no or only faint strips have developed 
yet; however, our focus here concerns later stages of the growth process. Then, fluctuations 
of the downstream interface position are rather small and very slow, so that its location can 
be determined very accurately. 

The evolution of the charge density profiles proceeds in three stages, as illustrated in Fig. 2 
for systems with different Ly. At early times, the downstream interface forms and equilibrates 
quite rapidly, reaching its steady state form at about 10^ MCS, independent of Ly. Near this 
interface, the charges are already perfectly segregated, with the charge profile saturating 
at ±1. Two fronts, one on each side, separate the ordered from the remaining disordered 
region. The larger systems now enter a second regime, during which both fronts travel slowly 
outwards, increasing the width of the ordered region. This process is also independent of Ly: 
the centers of the fronts move out as Ini, while their shape remains unchanged. Eventually, in 
the third regime, the two fronts merge, due to PBC, and the upstream interface equilibrates. 
The system has now reached steady state. For smaller systems (including the one shown in 
Fig. 1), the crossover to steady state sets in before well-developed fronts can form. 

As a quantitative probe, we define an order parameter, Q = ip'^{y,t) pc| ] which mea- 
sures the area, and hence the degree of order, under the (averaged, squared) charge profile. 
A central result of our work. Fig. 3 demonstrates that Q exhibits dynamic scaling in the 
logarithmic growth regime: Excellent data collapse is observed, for a range of E, Lx and Ly 
if EQ is plotted vs the scaling variable t = tE^ /L^. This confirms our picture of the ordering 
process: the motion and shape of the fronts is independent of Ly, and Q{t) oc Int, since each 
escape of the vacancy is an activated process (exponentially suppressed by E) . 

Two comments are in order. First, we should emphasize that Q differs from an order 

pa„™t„ „»d p„v.„„.l,, .a„e,y . (e. ( A E. »,) The ,att» ca„ be „e..„ed 

directly from the configurations, without the shifting procedure. Essentially, Q counts the total 
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Fig. 3 - Scaling plot of order parameter. The straight line, with slope 2.1 ± 0.1, is a fit to the data. 

number of ordered rows transverse to the external field. It is therefore sensitive to the presence 
of multiple strips which can emerge during the ordering process, especially in larger systems. 
As a result, the growth of Q exhibits a weak Lj,-dependence. In contrast, the shifting procedure 
averages out multiple strips, in favor of the largest (dominant) strip. Here, we focus on Q but 
note that the study of multiple strips remains an outstanding problem. Second, our scenario 
of the ordering process relies on having a well-developed downstream interface sandwiched 
between two fully charge-segregated regions. Our steady-state studies j20) show that this is 
the case provided ELy > 18. Thus, the system with Ly — 24, E — 0.8 sets a lower limit. 

Analytic approach. Finally, we turn to an analytic description, with the aim of deriving 
the appropriate scaling variables and the slope of the Int law. Our starting point is a set of 
mean- field equations for the (coarse-grained) local hole and charge densities. Due to particle 
number conservation, the continuum version of our model takes the form of two continuity 
equations which can be derived from a microscopic master equation. Their form has been well 
established |^,|l^,|2^. Proceeding directly to the equations for the profiles, we obtain 



Here, d denotes a spatial derivative in the drive direction, and terms of 0{(fP') have been 
neglected, since they reflect the presence of multiple vacancies. The continuum limit has been 
taken by letting the lattice constant vanish at finite E. The equations have to be supplemented 
with appropriate boundary conditions and the constraints on total mass and charge. 

Solving these coupled, nonlinear PDE's is of course highly nontrivial. However, the data 
suggest a few approximations which allow for considerable analytic progress. Since we wish 
to capture the logarithmic growth regime, we center the downstream interface at the origin 
{y — 0) and focus on the front moving in the positive y-direction (the other front being a 
mirror image, of course). By the onset of the second regime, this front is located well away 
from y ~ Q, since the charge density at the downstream interface has already equilibrated (cf. 
Fig. 2). Moreover, MC data show that the hole density profile is strongly localized at the origin 
and does not change significantly as time progresses. Therefore, we replace it in Eq. (||) by 
its steady state form (t>o{y) |^ which decays, to excellent approximation, as exp(— i?|?/|) 
on both sides. Here, c cx LxE~^ is a normalization, such that 0o(j/) is the probability of 
finding the vacancy at site {x, y) in the system. Introducing the scaling variables, z = Ey and 
T = c~^E^t, we arrive at a parameter-free equation, valid for y > 0: 



dt(t){y,t) 



d {d(j) + Etjytl;} 
d{(t)dil:-i:d(t)-E(t)} 



(1) 
(2) 



a,V(2, r) = exp(-z) [dli, - V + 1] • 



(3) 
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Clearly, initial and boundary conditions are needed. To study the motion of the front well 
before saturation occurs, we may allow < z < oo |2^. The disordered phase far ahead of the 
front is modelled by lim^^oo '>ij{z^ r) = for all r < oo. Further, the fully ordered phase well 
behind the front should be independent of i, since it is unaffected by further growth; thus, we 
demand that ?/'(z, t) ^ 1 for t ^ cxd and all < z < oo. 

Some comments are in order. First, since c oc Lx/ E, the characteristic scaling of MC time 
with E'^/Lx already emerges. Second, since Q = 2 ilj'^{y,t)dy = 2E~^ il;'^{z,T)dz, our 
mean field theory predicts data collapse if QE is plotted vs tE^/Lx, as borne out by Fig. 3. 

Since the moving front retains its shape, we may seek a solution in the form "0(2, t) = /(ui), 
with w = z — Zo{t). Clearly, Zo{t) describes the front position. Eq. (^) becomes 

- Zoe^"+^f' = /"-/ + 1, (4) 

where Zq = drZo and /' = dwf, with lim„^+oo f{w) = and lim„^_oo fiw) ~ 1. Also, /' 
differs notably from zero only in a (finite) neighborhood of w = 0. Multiplying both sides of 
(Q) by /', and integrating from — oo to +oo, we arrive at ioe^° dwe^ = i. Obviously, 

JJ*"^ dwe"^ f'^ = K is a positive, nonvanishing numerical constant. The time dependence of the 
front position (with simple initial condition Zo{0) — 0) now follows as Zo{t) = ln(l + t/2k). 
It is logarithmic, as expected. 

Returning to Eq. (|^), and using ioe^° — 1/2k, the shape of the front, subject to the 
specified boundary conditions, is given by f{w) = 2K,e~^ [1 — exp(--e"'/2K)]. The associated 
order parameter is now easily computed: QE = 21nr — 0.927... + 0(l/r), for large times 
T » K. Gratifyingly, the key features of the MC data, namely, the correct scaling variables, 
the logarithmic growth law and even its amplitude (i.e., the factor 2) are reproduced by our 
solution. Of course, a single additive constant is required to match the time scale of the 
simulations. 

Conclusions. In this letter, we presented MC and analytic results for domain growth in 
a simple driven lattice gas. A single vacancy rearranges positive and negative particles into 
charge-segregated strips transverse to the drive. Focusing on the largest strip, we find that it 
grows logarithmically with time, via field-suppressed excursions of the vacancy away from the 
center of the strip. The ordered domain is separated from the remaining disordered region by 
two well-defined moving fronts which retain their shape during the growth process. Eventually, 
the two fronts merge due to our BC's, and the system crosses over into steady state. During the 
logarithmic growth regime, we observe dynamic scaling of a suitably defined order parameter, 
Q, if QE is plotted vs tE^ /Lx- These findings are supported by an (approximate) solution of 
a set of mean-field equations for the charge and hole density profiles. To leading order in t, 
the scaling function QE grows as 2 In (tE^/Lx), in good agreement with the data. 

Several questions remain open. First, multiple-strip configurations form easily, especially 
in larger systems, and our focus on the dominant strip washes out any secondary ones, due 
to the shifting procedure. Unfortunately, a full dynamic scaling study of observables sensitive 
to multiple strips demands much greater computational effort. Whether similar coarsening 
behavior emerges remains to be explored. A second natural extension of our study would 
allow for multiple vacancies, i.e., a finite hole density. Under these conditions, particles 
first form "clouds" which coarsen and eventually merge into multiple strips transverse to the 
field. These strips then continue to coarsen until a single charge-segregated strip survives. A 
numerical solution of the ID mean-field equations p^ ] for the multiple-strip regime indicates 
that logarithmic growth persists, consistent with the arguments of Kafri, et.al. However, 
a detailed analysis, which identifies the relevant dynamic scaling variables and functions, and 
incorporates the crossover from clouds to strips, is still outstanding. 
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